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ABSTRACT 


The  Interactive  Signal  Synthesis  System  is  a  user- 
oriented  Fortran  program  designed  to  give  the  signal  analyst 
the  ability  to  generate  controlled  time  series  of  arbitrary 
length  and  complexity.  With  the  system,  one  can  not  only 
parameterize  basic  signal  functions  such  as  tone,  gate,  noise, 
polynomial,  and  arbitrary  data,  but  also  parameterize  the 
shifted  sum,  damping,  modulation,  and  filtering  of  each 
function. 


ADMINISTRATIVE  INFORMATION 

The  work  reported  here  was  performed  in  the  Computer  Science  and 
Information  Systems  Division  of  the  Computation,  Mathematics,  and  Logistics 
Department,  within  Program  Element  61153N,  Task  Area  SR0140301,  Task  15321, 
and  Work  Unit  1808-010,  the  Mathematical  Sciences  Research  Program  under  the 
sponsorship  of  NAVSEA  03R22. 


INTRODUCTION 

The  Interactive  Signal  Synthesis  System  is  a  user-oriented  Fortran 
program  designed  to  give  the  signal  analyst  the  ability  to  generate  controlled 
time  series  of  arbitrary  length  and  complexity.  With  the  system,  one  can  not 
only  parameterize  basic  signal  functions  such  as  tone,  gate,  noise, 
polynomial,  and  arbitrary  data,  but  also  parameterize  the  shifted  sum, 
damping,  modulation,  and  filtering  of  each  function. 

For  the  convenience  of  the  non-computer-oriented  analyst,  emphasis  has 
been  placed  on  structuring  the  interactive  user  interface  with  prompts, 
definitions,  and  optional  echoes.  However,  for  exact  formula  definitions  and 
discussions  of  side  effects,  the  user  must  refer  to  the  material  which  follows 
presently.  Output  signal  data  may  be  formatted  and  scaled  for  viewing  and 
manipulation  by  ISPARS  (Interactive  Signal  Pattern  Analysis  and  Recognition 
System).  Although  the  synthesis  system  has  been  coded  specifically  in 
Fortran  for  the  PDP-11/45,  it  contains  no  graphics  or  assembly  language  code; 
its  portability  therefore  is  strictly  a  function  of  its  machine-dependent 
input/output . 


*A  complete  listing  of  references  is  given  on  page  17. 


It  it  expected  that  the  signal  synthesis  system  will  be  of  use  in  generating 
test  data  for  validation  of  signal  analysis  programs,  in  establishing  an  analysis- 
by-synthesis  procedure  for  pattern  classification  applications,  and  in  determining 
standards  for  transient  waveform  definition. 

The  following  functional  description  of  the  signal  synthesis  system  outlines 
the  algorithmic  aspects  of  the  system.  A  separate  report  to  be  issued  as  an 
interactive  users'  guide  will  describe  the  decisions  required  of  a  user  in  executing 
the  program. 


FUNCTIONAL  DESCRIPTION 


INFORMATION  FLOW 

In  the  major  loop  of  the  Interactive  Signal  Synthesis  System,  a  sum  of  sampled, 
composite  functions  is  computed  as  follows: 

g(t)  -  J  f.(t)  (1) 

1 

where  the  time  t  is  related  to  the  sampling  period  Tg  by 

t  -  mTg  (m*0 ,  1 , . . . )  (2) 

The  reaultant  time  series  g(m)  is  placed  on  disk  in  INTFILE  in  integer  form  and/or 
on  SUMFILE  in  floating-point  form.  For  viewing  the  time  series  on  ISPARS,  an 
integer  input  file  is  required.  The  INTFILE  serves  this  purpose  and  is 
autosiatically  scaled  such  that  the  maximum  absolute  value  in  the  file  equals  8192. 

The  sum  g(t)  in  Equation  (1)  may,  in  fact,  be  represented  in  its  optionally 
filtered  form  g'(t)  as  follows: 

*,(t>  •  E,  <3> 

k-1 

where 

t'  -  t+((n+l)/2]T8-(k-l)T# 


t 


The  brackets  [,]  denote  "the  greatest  integer  in",  and  the  (k*l,...,n)  are  the 
weights  of  a  symmetric  FIR  (finite  impulse  response)  filter  to  be  described  more 
carefully  later.  Equation  (3)  is  then  the  convolution  of  a  filter  jw^}  with  a 
function  g(t)  transposed  in  time  so  that  the  center  k  ■  [Cn+l)/2]  of  the  symmetric 
filter  (which  corresponds  to  the  maximum  weight)  coincides  at  time  t  with  g(t).  The 
function  f^(t)  of  Equation  (1)  is  related  to  the  simpler  function  e^(t)  by  the 
shift-sum 


where 


R 

fi<t)  "  I  e.(t-(j-l)T  ) 

j-1 

R  (>0)  is  the  number  of  recurrences 
(>0)  is  the  recurrence  shift 
f^(t)  ■  0  for  t<0 


(4) 


By  judicious  use  of  the  damping  factor  to  be  associated  with  every  elementary 
function  to  be  described  below,  one  can  generate  a  repeating  sequence  of  disjoint  or 
overlapping  transient  signals  with  the  shift-sum  mechanism. 

The  shift-sum  f^(t)  of  Equation  (4)  may,  in  fact,  be  represented  in  its 
optionally  filtered  form  f|(t)  as  follows: 


fl(t)'"  1.  fi(t,,)  vp-k+l 
k-1 

where 

t"  -  t+([(p+l)/2j-(k-l))T8 


(5) 


Equation  (5)  represents  the  filtering  of  a  function  f^(t)  by  a  symmetric  FIR  filter 
with  weights  (k*l,...,p)  as  in  Equation  (3).  The  function  e^(t)  of  Equation  (4) 
may  then  be  formed  from  an  optional  modulation  M  of  the  parameters  of  the  elementary 
functions  d^(t).  Formally,  one  may  write 


e.(t)  -  e.(M,d.(t)) 


(6) 
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There  are  currently  five  basic  types  of  elementary  functions  d^(t): 

•  tone 

•  gate 

•  polynomial 

•  noise 

•  data 

each  of  which  will  be  defined  in  the  sections  to  follow.  The  modular  nature  of  the 
system  makes  it  easy  to  expand  the  list.  In  Equation  (6),  the  modulation  operator  M 
applied  to  d^(t)  is  implemented  by  the  use  of  some  d|(t)  selected  as  any  one  of  the 
five  listed  elementary  functions  from  which  d^(t)  is  itself  chosen.  Finally,  each 
d^(t),  d!(t)  is  subject  to  an  exponential  damping  which  is  included  in  the 
elementary  function  parameter  definition. 

In  terms  of  the  functions  just  briefly  described,  the  flow  of  information  and 
interplay  of  control  between  parameter  definition  and  time  series  synthesis  is 
illustrated  in  the  Figure  on  the  next  page.  In  the  Figure,  labelled  lines  signify 
information  and  unlabelled  lines  denote  control.  Thus,  the  execution  sequence  is: 
define  the  elementary  function  d  (and  its  damping),  its  modulation  type  M  (if  any), 
and  the  modulation  elementary  function  d'  (and  its  damping);  at  each  sample  m, 
evaluate  the  modulation  function  d1  and  apply  it  to  the  elementary  function  d  by  the 
operation  M(d);  shift-sum  and  optionally  define  and  apply  a  filter  \r.  The 
parameters  and  synthesized  signal  are  saved  on  their  respective  files.  Each 
iteration  will  sum  the  resultant  time  series  onto  the  INTFILE  and/or  the  SUMFILE 
with  past  syntheses  of  the  same  program  execution  and  permit  the  definition  and 
application  of  a  filter  w.  By  saving  the  parameter  file,  one  retains  the  ability, 
not  only  to  interactively  and  selectively  modify  a  parameter  set  previously  defined, 
but  also  to  resynthesize  a  signal  without  interaction.  The  latter  capability  is  a 
time-for-space  trade-off  that  permits  access  to  many  precisely  defined  time  series 
without  explicitly  storing  them. 


4 


PARAMETER  (PAR)  DEFINITION 
►ENTER 


F=> 


l 


T8,Tr,R 


PAR  function  (d) 
PAR(d) 

P 

PAR  modulation  (M) 


PAR(d,M) 

HT  PAR(d,M,d'),T  ,T  ,R 

PAR  function  (d1)  . . . -  . — 


PAR  filter  (v) 


PAR(v) 


I 


PAR(v) 

PAR  filter  (w) 


PAR(w) 


PAR  (w),li8t 


PARFILE 


ITERATE 


SIGNAL  SYNTHESIS 


Hr 

MOD  function  evaluate 


d’  (m) 


Hr 

MODULATE  d(m) 


e(m) 


Hr 

SHIFT- SUM 


f(m) 


J 

3  FILTER 


SUM 


J^g(m) 

3  FILTER 

r  »  a 


g'  (m) 


Hr 

INTFILE,  SUMFILE 


l 

EXIT 


Figure  -  Information  Flow  and  Control  Structure  of  Parameter  Definition 

and  Signal  Synthesis 
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PARAMETERS  OF  ELEMENTARY  FUNCTIONS 


In  this  section,  the  elementary  functions  listed  previously  are  defined,  with 
their  parameters  and  modulation  characteristics.  Each  description  is  self-contained 
and  therefore  some  redundancy  appears. 

Tone  Parameters 

The  exponentially  damped  tone  is  a  fundamental  function  from  which  any  single¬ 
valued,  periodic  time  series  of  engineering  interest  may  be  formed.  Nevertheless, 
certain  periodic  functions,  such  as  sawtooth  or  pulse  train,  which  require  a 
truncated  approximation  to  an  infinite  set  of  harmonics,  are  more  conveniently  and 
precisely  synthesized  by  more  direct  definitions,  such  as  the  gate  and  polynomial 
functions  to  be  described  later.  Moreover,  the  user  should  be  aware  that  a  multi¬ 
valued  function  such  as  a  true  gate  function  can  be  approximated  on  the  current 
system  only  by  a  function  whose  values  are  separated  by  at  least  one  sample.  The 
tone  function  is  defined  as  follows: 

d^t)  «  A  exp( -ALPHA  t)  sin(2lKFREQ  t+PHASE))  ♦  MEAN  (7) 


where 

A  is  the  amplitude 

ALPHA  (>0)  is  the  damping  constant  in  Hertz;  (ALPHA/2tt 
is  the  resonant  bandwidth) 

FREQ  is  the  frequency  in  Hertz;  (1/FREQ  is  the  period) 

PHASE  is  the  phase  as  a  fraction  of  2t r;  i.e.,  if  P  is 
the  phase  in  radians,  then  PHASE  ■  P/27 r 
MEAN  is  the  mean 

One  and  only  one  of  the  parameters  amplitude,  frequency,  or  phase  is  subject  to 
modulation,  which  has  the  form 


Pt(t)  -  p[d[(t) 


(8) 
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where 

Pj(t)  is  one  of:  A(t),  FREQ(t)/t,  PHASE(t) 
pj  is  the  value  of  p^  chosen  when 
d^(t)  was  defined 

d!(t)  is  one  of  the  five  elementary  functions 

The  option  of  simultaneous  modulation  modes  is  not  implemented  but  is  under 
consideration.  The  user  should  note  that  the  frequency  modulation  defined  by 
Equation  (8)  is  spectral  frequency  modulation,  not  the  more  familiar  instantaneous 
frequency  modulation.  In  order  to  produce  the  latter  type  of  modulation,  one  need; 
to  implement  the  following  relation: 


t 

FREQ(t)  «  FRfQ  ( 

0 

so  that  the  phase  derivative  -  the  instantaneous  frequency  -  is  the  desired  function 
d!(t).  One  can  implement  Equation  (9)  by  synthesizing  d!(t),  passing  it  through  a 
low-pass  filter  as  described  later  in  the  section  on  filter  parameters,  and 
frequency  modulating  d^(t)  by  the  resultant  floating-point  SUMFILE  (plus  a  CONSTANT) 
regarded  as  a  function  d^(t),  as  explained  later  in  the  section  on  the  data  function 
parameters.  If  d^(t)  is  itself  chosen  to  be  a  tone  d|(t),  then  instantaneous 
frequency  modulation  is  implemented  simply  as  phase  modulation  by  delaying  the 
modulating  tone  by  w/2  and  multiplying  by  the  reciprocal  of  2ttFRQ' ,the  radial 
frequency  of  the  modulating  tone. 

Gate  Parameters 

As  a  means  for  introducing  logical  gating  sequences  and  gate  train  synthesis, 
the  gate  function  is  now  defined. 


/ 


d^(t)  dt  +  CONSTANT) 


(9/ 
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T 


d2(t) 


(OFFSET  A  exp( -ALPHA  (t-PHASE  Tp) ) , 
t-PHASE  Tp  <  Tg  (mod  Tp) 


^  OFFSET ,  t-PHASE  Tp  >  TQ  (mod  Tp) 


(10) 


where 

OFFSET  is  the  constant  offset  (i.e.,  the  range  mean,  not 
the  function  mean) 

A  is  the  offset  to  peak  amplitude 
ALPHA  (>0)  is  the  damping  constant  in  Hertz 
PHASE  is  non-negative  phase  (i.e.,  delay)  as  a  fraction  of 
the  period  Tp  ;  i.e.,  if  D  is  the  delay  in  seconds 
(mod  Tp) ,  then  PHASE  -  D/Tp 
Tp  is  the  period;  gate  frequency  is  FRQ  ■  1/Tp 
Tg  is  the  gate  in  seconds  (i.e.,  the  on-period  duration,  not 
greater  than  Tp) ;  duty  cycle  is  Tg/(TG+Tp) 

One  and  only  one  of  the  parameters  amplitude,  phase,  period,  or  gate  is  subject 
to  modulation,  which  takes  the  form 

P2(t)  *  p£d[(t)  (11) 

where 

P2(t)  is  one  of:  A(t),  PHASE(t),  Tp(t),  Tg(t) 

P2  is  the  constant  value  of  p2(t)  chosen  when  d2(t) 
was  specified 

d^(t)  is  one  of  the  five  elementary  functions 

The  option  of  simultaneous  modulation  modes  is  not  implemented  but  is  under 
consideration.  Note  that  if  Tg*0,  a  single  sample  is  defined  as  the  on-period,  thus 
producing  an  impulse  train. 
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Polynomial  Parameters 

A  generalized  form  of  the  gate  function  for  shaped  pulse  trains  may  be 
synthesized  by  a  polynomial  function  of  order  q  as  follows: 


d3(t) 


/  OFFSET  +  A  exp(-ALPHA(t -PHASE  T  )  V  a.  (t-PHASE  T_) 1 , 

I  F  i-0  1  F 


I 


OFFSET, 


t-PHASE  T„  <  T_ 

r  —  u 


t-PHASE  T  >  T_ 
F  G 


(mod  Tp) 


(mod  Tf) 


(12) 


where 


OFFSET  is  the  constant  offset 
q  is  the  polynomial  order 

a^  (i*l,...,q)  are  the  polynomial  coefficients 
A  is  the  offset  to  peak  amplitude 
ALPHA  (>0)  is  the  damping  constant  in  Hertz 
PHASE  is  non-negative  phase  (i.e.  delay)  expressed  as  a 
fraction  of  the  period  T  ;  i.e.  if  D  is  the  delay 


in  seconds  (mod  Tj.) ,  then  PHASE 


D/Tt 


Tp  is  the  period  in  seconds;  frequency  is  FRQ  *  1/Tp 
T  is  the  gate  in  seconds  (i.e.  the  on-period  duration,  not 
greater  than  Tf) ;  duty  cycle  is  TG/(Tg+TF) 

One  and  only  one  of  the  parameters  amplitude,  phase,  period,  or  gate  is  subject 
to  modulation,  which  has  the  following  form: 


P3(t) 


p^d 


(13) 


where 

P3(t)  is  one  of:  A(t),  PHASE(t),  TF(t),  Tg(t) 

P3  is  the  constant  value  of  p^(t)  chosen  when  d^(t) 
was  specified; 

d|(t)  is  one  of  the  five  elementary  functions 
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When  Equation  (12)  ia  applied  to  to  the  generation  of  a  polynomial  of  specific 

order,  care  must  be  taken  to  choose  the  a^  properly.  Some  examples  follow.  In  each 

example,  the  function  will  start  at  a  level  X  at  t„,  and  turn  off  at  t_  by 

on  U  T 

switching  to  the  le„vel  Xo^.  Thus,  from  Equation  (12),  one  can  immediately  write 


*0  ■  (Xo«-Xoff)/2 

( 14) 

OFFSET  ■  (X.  .X  ,,)/! 

on  ocr 

(15) 

Example  1.  Ramp  sequence.  -  If  a  first  order  polynomial  is  chosen  as  the  basic 
pulse  shape,  a  gated  ramp  sequence  results.  The  ramp  is  defined  as  rising  or 
falling  to  a  target  value  Xtflr  in  a  duration  time  T.  .  The  equation  for  the  first 
order  coefficient  is  then  as  follows: 


a,  •  (X_  -X  )/T, 

1  tar  on  dur 


(16) 


Example  2.  Triangle  sequence.  -  A  gated  triangular  wave  results  from  the 

synthesis  of  the  sum  of  two  ramp  sequences.  If  the  peak  has  a  value  Xtar  which  is 

reached  in  duration  time  T^^  from  the  origin,  then  the  value  for  a^  is  the  same  as 

in  Equation  (16).  If  the  triangle  terminates  at  a  terminal  value  Xfcer  at  time  sTdur 

from  the  origin,  then  the  second  ramp  begins  at  time  T.  at  level  zero  and  has  the 

dur 

slope 


<X.«r-'Xt.r*<,-1)lton)«*-»Tdur 


07) 


Example  3.  Parabolic  sequence.  -  A  third  order  polynomial  chosen  as  the  basic 
pulse  shape  will  produce  a  gated  parabolic  sequence.  If  the  parabola  rises  or  falls 
to  successive  values  Xmed,  Xtar  at  time  intervals  T^,  Tj,  respectively,  the 
equations  for  the  two  higher  order  coefficients  are  as  follows:  (Note  that  the  total 
gate  duration  Tdur  is  the  sum  of  the  intervals  T^ ,  Tj  just  mentioned.) 
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*1  ■  ^Tl^Xon-Xtar)+^dur^«ed”*on^^T1^2^Tl+^2^ 


a.  -  (T.Xt  +T-X  -T.  X  .)/T,T,(T,+T») 
2  1  tar  2  on  dur  med  121  2 


(18) 

(19) 


If  the  gated  parabola  sequence  were  being  used  to  frequency  modulate  a  tone,  the 
coefficient  a^  would  have  to  be  adjusted  for  the  fact  that  the  instantaneous 
frequency  would  be  computed  by  differentiating  the  polynomial,  thus  doubling  the 
contribution  of  the  a^  term. 


Noise  Parameters 

A  representation  of  random  behavior  may  be  introduced  in  the  synthesized  signal 
by  drawing  samples  from  either  uniform  or  Gaussian  distributions. 

In  the  uniform  case,  a  random  variable  RAN  uniform  on  (0,1>  is  derived  from  a 
FORTRAN  library  function  that  requires  two  user-supplied  integers  to  initiate  the 
generation.  The  uniform  noise  function  with  mean  MEAN  and  standard  deviation  A 
SIGMA  ia  then  defined  as  follows: 


d4u(t)  -  A  SIGMA  ^12  exp( -ALPHA  t)  (RAN(t)-.5)  ♦  MEAN  (20) 


where 

A  is  the  amplitude  gain  '*> 

ALPHA  (>0)  is  the  damping  constant  in  Hertz 
RAN(t)  is  a  uniform  random  variable  on  (0,1>  computed 
recursively  from  RAN(t-l) 

SIGMA  is  the  standard  deviation  for  A* 1 
MEAN  is  the  mean 

In  the  Gaussian  case,  a  univariate  distribution  may  be  derived  from  two 

.  2 

independent  uniform  distributions,  as  explained  in  a  paper  by  Franklin.  The 
appropriate  time  series  is  divided  into  even  and  odd  components  as  follows: 
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(21) 


d4a<2t)  -  A  SIGMA  exp (-ALPHA  2t) 

^-2  ln(RANA(t))  cos(RANB(t)) 

♦  MEAN 

d4jt(2t+l)  -  A  SIGMA  exp(-ALPHA(2t+l) ) 

^TT^rRANATtTy  sin(RANB(t )) 

♦  MEAN  (22) 


where 

A,  SIGMA,  ALPHA,  MEAN  are  as  above  for  uniform  noise; 

RANA(t),  RANB(t)  are  independent  random  variables  on  (0,1>, 
computed  recursively  from  RANA(t-l) ,  RANB(t-l), 
respectively. 

Only  amplitude  modulation  is  available  for  the  noise  functions,  in  the  form 

A(t )  -  A'df(t)  (23) 


where 

A'  is  the  defined  amplitude  of  d4(t) 

d!(t)  is  one  of  the  five  elementary  functions 

Data  Parameters 

The  final  elementary  function  is  in  a  sense  arbitrary.  Samples  generated 
either  empirically  or  by  previous  syntheses  may  be  read  from  any  integer  or 
floating-point  user-defined  data  file.  This  file  must  have  been  created  using  the 
PDP  11/45  SYSLIB  function  IWRITW.  The  form  of  the  data  function  is  as  follows: 

dj(t)  -  A  DATA(t)  +  OFFSET  (24) 


where 

A  is  the  amplitude  gain 
DATA(t)  is  the  data 
OFFSET  is  the  offset 


12 


Only  amplitude  modulation  ia  available  for  the  data  function,  in  the  form 


A(t)  -  Ad[(t) 


(25) 


where 

A  ia  the  defined  amplitude  gain  of  dj(t) 
d!(t)  is  one  of  the  five  elementary  functions 

FILTER  DESIGN 

As  indicated  in  the  Figure,  the  signal  being  synthesized  stay  be  filtered  after 

shift-sum  or  after  the  iterate  sum.  The  total  filter  may  be  a  composite  of  up  to  30 

lowpass,  bandpass,  bandstop,  or  highpass  filters,  each  of  whose  weight  vectors  w^  is 

designed  using  canned  software  written  by  Rabiner,  McGonegal,  and  Paul.  This 

3 

software  has  been  published  by  the  IEEE  in  a  more  comprehensive  filter  design 
package  than  that  offered  here.  The  composite  filter  vector  is  formed  by  a 
normalized  sum  of  the  component  filter  weights  w„  against  gains  g^,  as  follows: 


w 


n  n 

£  8 A.  1  £  *i 

l- l  1-1 


(26) 


where  the  filter  length  n  is  an  integer  lying  between  3  and  896.  For  highpass  or 
bandpass  filters,  n  must  be  odd.  The  component  filters  are  symmetric,  linear  phase 
FIR  digital  filters  designed  on  the  basis  of  cut-off  parameters  and  against  the 
following  windows: 

•  rectangular 

•  triangular 

•  Hanning 


•  generalized  Homing 

•  Hanning 

•  Chebyshev 

•  Kaiser-Bessel 

r 

3 

The  exact  forms  of  the  windows  may  be  found  in  the  IEEE  document  mentioned  earlier, 
page  5.2-2.  Here,  however,  the  parameters,  if  any,  associated  with  each  window  and 
their  effects  are  only  briefly  summarized.  Notwithstanding  the  theory,  the  user 
will  find  it  necessary  to  experiment  with  window  selection  and  parameter  values  to 
meet  his  needs.  Note  that  specification  of  a  desired  gain  parameter  and  appropriate 
frequency  cutoffs  is  common  to  all  the  filter  designs. 

The  rectangular,  triangular,  Hamming,  and  Hanning  windows  have  no  user-defined 
parameters  other  than  gain  and  cutoffs. 

The  generalized  Hamming  window  has  one  parameter,  ALPHA,  ranging  over  <0,1>. 

The  following  table  illustrates  the  effect  of  varying  ALPHA: 


TABLE  1  -  GENERALIZED  HAMMING  WINDOW  PARAMETER  VARIATION 


ALPHA  - 

WINDOW  - 

1.0 

rectangular 

0.54 

Hanning 

0.5 

Hanning 

0.0 

cos(2rri/n)  for  -n  £  4i  n 

For  the  Chebyshev  window,  in  addition  to  specifying  the  frequency  cutoff  and 
desired  passband  gain,  one  must  specify  two  of  three  parameters  from  which  the  third 
is  calculated.  The  parameters  are:  filter  length  n,  transition  width  DF,  and  ripple 
DP.  The  transition  width  is  the  length  of  a  frequency  band  between  the  passband  and 
stopband.  These  bands  are  defined  as  follows.  Let  G  denote  the  desired  passband 
gain.  Then  the  frequency  for  which  the  computed  gain  first  falls  below  G-DP  is 
called  the  passband  cutoff.  Similarly,  the  frequency  for  which  the  gain  first  falls 
to  DP  ia  called  the  stopband  cutoff.  The  difference  between  the  passband  and 
stopband  cutoffs  is  called  the  transition  width  DF.  It  is  the  passband  cutoff  which 
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will  have  been  specified  as  the  desired  frequency  cutoff  by  the  user.  The 
passband  ripple  is  defined  to  be  the  maximum  deviation  in  decibels  above  the  \ 

desired  gain  in  the  passband.  Similarly,  the  stopband  ripple  is  defined  to  be 
the  deviation  from  zero  gain  in  the  stopband.  In  the  present  design, 

both  passband  and  stopband  ripples  are  defined  equal  to  DP.  For  given  filter 

length  n,  DF  and  DP  vary  inversely. 

The  Kaiser-Bessel  window  has  one  parameter  ATT  which  is  the  desired 

minimus  stopband  attenuation  in  decibels. 
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